Quantum current noise from a Born-Markov master equation 
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Quantum coherent oscillations in the electric charge passing through a mesoscopic conductor can 
give rise to a current noise spectrum which is strongly asymmetric in frequency. The asymmetry 
reveals the fundamental difference between quantum and classical fluctuations in the current. We 
show how the quantum current noise can be obtained starting from a Born-Markov master equation, 
an approach which is applicable to a wide class of systems. Our method enables us to analyse the 
rich behavior of the current noise associated with the double Josephson quasiparticle resonance in 
a superconducting single-electron transistor (SSET). The asymmetric part of the noise is found to 
be strongly dependent on the choice of operating point for the SSET and can be either positive or 
negative. Our results are in good agreement with recent measurements. 

PACS numbers: 72.70.+m, 42.50.Lc, 85.25.-j 



Introduction. Electrical circuits have long provided a 
test-bed for studying fluctuations and noise— Recent 
attention has focussed on the differences between quan- 
tum and classical noise far from equilibrium^. Classi- 
cally, the correlation function describing the average of 
the product of currents measured at two times is real and 
symmetric with respect to time. In contrast, in quantum 
mechanics current is described by an operator which need 
not commute with itself at different times, leading to 
quantum current correlation functions that are complex 
and asymmetric. The asymmetry in quantum correlation 
functions ultimately arises because the probabilities of a 
quantum system emitting or absorbing a given amount of 
energy are not in general the same-4 The spectral den- 
sity of the quantum current fluctuations, S7(w), is real 
but asymmetric in frequency and can be measured by 
coupling the circuit to a mesoscopic detector and measur- 
ing the energy absorbed and emitted— ~—. Asymmetry in 
Si(cj) has been observed in superconducting tunnel junc- 
tions— quantum point contacts—^, a carbon nanotubc 
quantum dot in the Kondo regime— and superconduct- 
ing single-electron transistors (SSETs)— ~— . 

A common approach to describing transport in meso- 
scopic devices, especially when Coulomb blockade effects 
are important, is to derive a master equation for the den- 
sity matrix of the system by starting from the Hamil- 
tonian evolution of the system coupled to leads which 
act as an environment. For many devices, such as quan- 
tum dots 14 d 5 or single electron transistors— ~—, it is then 
possible to use the Born-Markov approximations to de- 
scribe the system using an equation of the form, 



p = £p, 



(1) 



where p is the reduced density operator obtained by trac- 
ing over degrees of freedom associated with the leads and 
jC is a super-operator—. 

The calculation of current correlation functions is more 
complicated than for other quantities, such as the charge 



accumulated in a dot or a spin projection, because of 
the difficulty of defining a current operator. Current in- 
volves keeping track of particles in the leads which do not 
necessarily form part of the system in a Born-Markov de- 
scription. For this reason many calculations have consid- 
ered just the symmetric (classical) current noise 18 i 24 ~ — , 
or used equations of motion for auxiliary operators to 
obtain symmetrized correlation functions—. Those cal- 
culations that do address the quantum noise have con- 
centrated on effects which are important on very short 
time-scales where the Markov approximation fails— ~—. 

In this Rapid Communication we show how the quan- 
tum current correlation functions can be obtained from a 
Born-Markov description. This approach provides a way 
of calculating the quantum current noise arising from co- 
herent oscillations in the charge passing through a con- 
ductor and can be applied to a range of systems, includ- 
ing the SSETs probed in recent experiments— ~—. Us- 
ing our method, we provide a theoretical analysis of the 
quantum current noise spectrum near the double Joseph- 
son quasiparticle (DJQP) resonance. Our results reveal 
a complex behavior with regions of positive and nega- 
tive asymmetry in the noise, depending on the precise 
choice of SSET operating point, in accord with measure- 
ments—. 

General Formalism. Wc would like to calculate the 
fluctuations of the current between a given lead and 
the system itself for a device whose dynamics is de- 
scribed by Eq. ([1]). A convenient basis for the Hilbert 
space is {\N,n)}, where N\N,n) = N\N,n), with N 
an operator counting the number of charges in the lead, 
and n labeling the other quantum numbers needed to 
identify a given state. In this basis the density ma- 
trix has elements: p a {N) = (N,n\ p\N + q,n') with 
a = {n, n';q}. Non-diagonal elements of p with q 
can be non-zero when transport is coherent, for instance 
in the presence of a Josephson junction. The equation 
of motion for p a (N) reads p a (N) = J2p, p £apPp( N + V), 
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where we used the fact that L cannot depend on N ex- 
plicitly. This equation is solved by Fourier transform- 
ing in N. Defining p a {x) = Hn ^ lxN Pa(N) one finds: 
P(.X>t) = e V(Xi 0), where p is a vector with compo- 
nents p a and C the matrix £ Q( g = J2 P s~ lxp C 



(p) 

a/3" 



Given p one can obtain the (time-dependent) average 
of any operator and, in particular, of N: 



(N(t)) = Tr{N(t)p} = w 



dix 



(2) 



where we introduced the transpose of the vector w, whose 
matrix elements w a arc 1 when a — {n, n; 0}, that selects 
the diagonal elements of p. The particle current is ob- 
tained by taking a time derivative (I(t)) = —d(N(t))/dt: 



(Ht)) 



w 



de c{x)t 
dix 



£'(x)e 



PO 



. (3) 

Here po is the initial density matrix, and a prime indi- 
cates a derivative with respect to %x- The first term of 
Eq. ([3]) vanishes since w t C(0) = 0. This relation fol- 
lows from the conservation of probability w p(0,t) = 1 
and the equation of motion for p(x,t). For the second 
term, the hypothesis that the process is stationary im- 



plies that lim t 



a £(0)t 



Pq = Pst f° r an y initial vector pq. 



Note, however, that in general p a (N,t) does not reach a 
stationary value, only p a (N,t) does. The resulting 
expression for the stationary current is then: 



(I) = -w'C'Mps 



(4) 



We now consider the current correlation function, 
Siituta) = (I(h)I(h)) - {1(h)) (1(h)). The first term 
can be written as follows: 

</(<i)/(*2)> = ^-^Tt{N(h)N(t 2 )p(t = 0)} • (5) 

We now define for h > h f(h,h) = (N(h)N(t 2 )). 
Since (N(t 2 )N(h)) = f(h,t 2 )*, Si(h,t 2 ) = d tl d t2 {6(h- 

h)f(h,h) + 0(h - h)f(h,h)*} - (1(h)) (Hh)), with 

9(t) the Heaviside step function. We can now use the 
quantum regression theorem— to obtain the correlator 
f(h,t 2 ) for h > h- 

f(h,h) = Tr ^Nc t{tl - t2) Nc £t2 p(t = 0) j , (6) 

where N is the super-operator that for any op- 
erator, A, has the property: {N, n\ NA \M, n') = 
N (N, n\ A \M, n'). Finally, we express Eq. (|6]) as a scalar 
product, 



f(h,h) = w l 



d 

dix 



a £( X )(ti-*2) 



d_ 

dix 



Po ■ 



(7) 



Performing the time derivatives we find that Si(h-,h) 
only depends on t = h — h for h — > oo. For t > the 
correlation function reads: 



S/(t) = SWBeiw'C'pst] + w t C'e ct C' Pst - (I) 



(8) 



where the x = argument is omitted for brevity (the 
corresponding result for t < is obtained by complex 
conjugation). The time Fourier transform of Eq. <j8j> is 
obtained by introducing the right and left eigenvectors of 
£(0): £(0)iv = \ v v v and ui*£(0) = A„u£, with = 
In particular, Ao = with wq = w, vq = p st , and 



ReA„ < for v > 0. Wc find: 



Si(uj) = Re 



w t C'v L , wjCps 



(9) 



Equation ^ is the central result of this paper: given a 
form for C(x), it can be used to compute the correspond- 
ing current noise spectrum. 

At first sight, Eq. ((9]) looks rather similar to expres- 
sions obtained by calculating probabilities for given num- 
bers of chargeSjJV, to have passed into the leads, such as 
in Refs.Hfjl andll^- However, there is a crucial difference: 
the calculation in Refs.[26| and[13 assumes that there is no 
coherence between states with different TV. Our choice of 
basis for p a (N) allows us to express as a simple deriva- 
tive with respect to ix, leading to Eq. © which does take 
fully into account coherence in N . Including such coher- 
ences is inherently problematic for methods based on the 
coun ting statistics of charges in the leads as discussed in 
Refs.[Tjand[3l 

Quantum Current Noise in SSETs. We now apply our 
method to the concrete example of the DJQP resonance 
that occurs in a SSETi2>2£3£. A SSET consists of a 
superconducting island coupled by Josephson junctions 
to superconducting leads; a voltage applied to a gate is 
used to tune the island potential. When a bias voltage 
is applied across the device, a combination of coherent 
Cooper-pair oscillations and quasiparticle tunneling can 
give rise to a stationary current. The DJQP resonance 
occurs for voltages where both resonant Cooper-pair tun- 
neling and quasiparticle tunneling occur at both junc- 
tions— The classical current noise for the DJQP 
has been calculated— (though only in certain limits). A 
recent experiment that probed the asymmetry in the cur- 
rent noise— provides motivation for us to study the full 
quantum problem. 

The Hamiltonian of the SSET is the sum of two terms: 
H = He + Hj. The charging part is 

H c = ^[E C (n - n g ) 2 + (N + n/2)eV] \N, n) (N, n\ , 

N.n 

(10) 

where V is the bias voltage, Ec is the island charging 
energy and n g is the number of island charges induced 
by the gate voltage (the SSET is assumed to be symmet- 
ric). The states {\N,n)} form a complete basis with N 
and n the number of electrons in the left lead and the 
island, respectively. The number of electrons in the right 
lead is — N — n plus a constant (which we set to zero). 
The Josephson part of the Hamiltonian relevant for our 
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FIG. 1: Schematic representation of the sequence of states 
involved in the DJQP cycle. The numbers give the electron 
occupation of the island and of each lead. 



problem is 

Hj = -jY, (W, °> ( N ~ 2, 2| + |JV, 1) (N, -1| + h.c.) , 

N 

(11) 

where J = Ej/2 with Ej the junction Josephson en- 
ergy. A DJQP resonance occurs for voltages such 
that n g = 1/2 and eV = 2Ec- the pairs of states 
{]iV, 0) , |iV — 2,2)} and {\N, 1) , \N, -1)} are resonant, 



while quasiparticle decays are possible (provided Ec > 
2A/3) between the states \N,2) -> \N,1) through the 
right junction and between \N, —1) — > N — 1, 0) through 
the left junction. (Other quasiparticle decays are blocked 
for Ec < 2A.) The sequence of transitions and corre- 
sponding changes in N and n are shown in Fig. [TJ 

Assuming that the normal resistance of the junc- 
tions, Rj, satisfies Rj ;§> h/e 2 , the dynamics is de- 
scribed by Eq. ([1]) with Cp = (C co h + Cdec)p, where 
C-cohP = —i[H,p] gives the coherent evolution. The 
term C decP = £ s T s [\s') (s\ p \s) (s'\ - (l/2){|s) (s\,p}}, 
describes dissipative quasiparticle tunneling de- 
cay process between states \s) and \s'). The states 
s) — > \s') concerned are \N, — 1) — > | TV — 1 , 0) with 
T s = T L and \N,2) -> \N,1) with r s = T R . For 
given TV, a set of 8 matrix elements of p is sufficient 
to describe the system 18 i 35 . Using the notation intro- 
duced above we introduce the vector: p = {po,0;0, P2,2;0> 

P0,2;-2, P2,0;2, P-1,-1;0, Pl,l;0, P-1,1;0, Pl,-l;o}- Writ- 
ten like this, the evolution equation for the first el- 
ement is p ,0;o( N ) = M [P2,0;-2(N - 2) - p ,2;2(N)] + 

rLp_i._i : o(-/V + 1). Fourier transforming with respect 
to N gives the first row of the matrix C(x)'- 
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(12) 



where the dctunings from the resonances for the left and 
right junction arc defined as Sl = ^Ec(n g — 1) + eV 
and S R = AEcn g — eV. Now applying Eq. ([§]) with w = 
{1,1,0,0,1,1,0,0}, numerical diagonalization gives the 
frequency dependence of the current noise through the 
left junction, Sl(w) for a given set of system parameters. 

Our quantum calculation reveals the asymmetry in the 
current-noise spectrum as well as giving insights into the 
high frequency behavior. Figure [2] shows the Fano fac- 
tor for the left junction, Fl{uj) = Sl{u)/(I). Defin- 
ing r as the quasiparticle decay rate at the center of 
the DJQP resonance (where Tl = T R ), we can dis- 
tinguish two different behaviors corresponding to weak 
(T/Ej <C 1) and strong (T/Ej 3> 1) quasiparticle tun- 
neling. For T/Ej -C 1 [sec Fig. [2fa)], coherent Cooper- 
pair oscillations lead to strong peaks in the spectrum at 
lu ~ iEj. At linear order in <5l(_r) the peaks have heights 
F L (±Ej) = 16Ej(Ej T 2S L )/(9T 2 ) and width oc T. For 
Sl < (5l > 0) the positive (negative) frequency part of 



Fl(uj) is enhanced since resonant oscillation in the SSET 
involves absorption (emission) of energy. The value of 
Sr influences the magnitude of Fl(uj), but not its asym- 
metry. This is because Sr affects the flow of current, 
but not the relative probabilities of energy absorption 
and emission at the left junction. For T/Ej ^ 1 [Fig. 
[2^b)] , there are no peaks in Fl{ui), but only a dip around 
w = with F L = 3/2 at \u>\ < E 2 /T and F L = 5/3 at 
E 2 j/T < |w| < r. 

For arbitrary T/Ej and at linear order in #£,(m, ui one 
finds again that the asymmetry is controlled uniquely 
by 5 L : F™ ym (uj) = (2/3)S l uj/(T 2 + 4Ej), where we use 
the notation asym (w) = 0(w) - <j>{-u). Note that the 
asymmetry has a purely quantum nature, in contrast to 
the symmetric part of the low frequency spectrum, which 
could have been obtained using the methods of Rcf. 

For frequencies to S> Ej,T we find Fl(lj) — > 1/3 in- 
dependent of all system parameters (though the Born- 
Markov approach breaks down eventually in the limit of 
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FIG. 2: (Color online) Frequency dependent Fano factor, 
Fl{uj). (a) Weak quasiparticle tunneling regime Y / E.j = 0.2: 
Sl = (black solid line), Sl — Y (blue dashed line) and 
Sl = — r (red dotted line), (b) Strong quasiparticle tun- 
neling regime Y/Ej — 50: Sl = (black solid line) and 
Sl = 0.2F (blue dashed line). In all cases Sr = and we 
set r_L = Tr = r neglecting the voltage dependence of the 
rates. 



very high frequencies L) 3> Ec)- This sub-Poissonian 
noise arises because Cooper pairs contribute to the cur- 
rent, but not to the high frequency noise. The asymmetry 
at high frequency for small $l(r) appears with the term 

M = -(4/3)(r 2 + 4^)<w^. 

The noise measured by coupling a detector to one of 
the leads is given by a combination of the particle and 
displacement currents— 



^asym 



Sj(w) = [S L (u) + S r (lj)}/2 - lj 2 S q (uj)/4, 



(13) 



where Sq is the charge noise spectrum (studied in 
Ref.l20h and we assume a gate capacitance much less than 
those of the junctions. The current noise at the right 
junction, Sr(co), is obtained using the same technique 
as for Sl(w) but using states which track the charge in 
the right lead. Using a similar method, we obtain an 
expression for Sq(u>) analogous to Eq. 



Sq(u) 



w t hv l/ wlnp st \ 
iuj + \ u J ' 



here h is a matrix corresponding to the super-operator 
representation of the island charge operator. 

In the experiments of Xue et al.— an electrical res- 
onator is used to probe asymmetry in the SSET cur- 
rent noise. Within the linear-response regime—, the 
SSET leads to damping of the resonator at the rate 
7(wo) = A 2 Sj Sym (u>o) , where A is the strength of the 
SSET-rcsonator coupling and ujq the resonator frequency 

Using the parameters in Ref.[l3l and Fermi's golden rule 
for calculating the quasiparticle rates— we obtain the 
contributions to 7 from each of the terms in Eq. (fT3")) . The 
particle current contributions, r )L(R) = A 2 S L ^^(u>o), 
shown in Figs. Eta) andE^b), are (almost—) antisym- 
metric about lines where the corresponding junction has 
a Cooper-pair resonance. Regions of positive (nega- 
tive) damping arise when a resonance is detuned so 
energy is on average absorbed (emitted) from the res- 
onator by the SSET. The charge noise contribution, 
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FIG. 3: (Color online) Contributions to the damping (in 
MHz) from particle currents in (a) left, 7l, and (b) right, "/r, 
junctions; (c) contribution from the charge noise, 7q and (d) 
total damping, 7. Dashed lines indicate the two Cooper-pair 
resonances which overlap at the DJQP resonance. Parame- 
ters used are: E c = 237 peV, Ej = 51 ^eV, A = 190 neV, 
Rj = 27kft, luo/2-k = 1.04 GHz and A = 2.28 x 10 18 C" 1 -^, 
corresponding to Y/Ej ~ 1.98 at resonance. 
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(u>o), [Fig. E^c)] has a different sym- 
metry as both Cooper-pair resonances affect the island 
charge. The overall damping, 7 = (7^ + Jr)/2 — 7q/4, 
shown in Fig. OJd), is dominated by 7^ and 7^; the in- 
fluence of 7q is weak because the frequency scale for the 
SSET is set by Ej, which is much larger than ujq. 

A simple comparison can be made with Ref.[l3| by com- 
puting the maximum and minimum values of 7, which 
occur for n g = 0.5 and bias voltages below and above 
the center of the DJQP resonance [see Fig. Eld)]. We 
obtain maximum and minimum damping rates with the 
same magnitude, 475 MHz, but opposite sign, in accord 
with the symmetry of the problem. Measured maxi- 
(14) mum and minimum damping rates— were 550 MHz 
and « —35 MHz respectively. Our calculation fits with 
the experiment on the low bias side, though agreement 
is less good on the high bias side. The difference is prob- 
ably due to the low resistance junctions (Rj = 27kf2) 
used 13 ' 34 which allow higher-order processes beyond the 
DJQP whose contribution to the current (and hence to 
the damping) increases with the bias voltage. 

Conclusions. We have shown that quantum current 
noise in a mesoscopic conductor can be calculated using 
a Born-Markov master equation description. The theory 
presented allowed us to find the asymmetry of the quan- 
tum current noise at the DJQP resonance in SSETs and 
to confirm the interpretation of a recent experiment that 
measured this asymmetry by detection of emission and 
absorption of energy The method we derived here has a 
wide scope of applicability. It could, for example, be ap- 
plied to Cooper pairs resonances in the SSET, for which 
theoretical predictions have not yet been made though 
the quantum current noise was measured recently—. 
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